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ABSTRACT 



We present some theoretical predictions concerning the amphtude and mag- 
netic sensitivity of the hnear polarization signals produced by scattering pro- 
Q> ■ cesses in the hydrogen Lya line of the solar transition region. To this end, we 

C/^ ■ have calculated the atomic level polarization (population imbalances and quan- 

1^ ■ tum coherences) induced by anisotropic radiation pumping in semi-empirical and 

I . hydrodynamical models of the solar atmosphere, taking into account radiative 

^ . transfer and the Hanle effect caused by the presence of organized and random 

^ . magnetic fields. The line-center amplitudes of the emergent linear polarization 

signals are found to vary typically between 0.1% and 1%, depending on the scat- 
tering geometry and the strength and orientation of the magnetic field. The 
results shown here encourage the development of UV polarimeters for sounding 
OO . rockets and space telescopes with the aim of opening up a diagnostic window for 

—I. . magnetic field measurements in the upper chromosphere and transition region of 

[-^ . the Sun. 
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1. Introduction 

One of the big cliallenges of 21st century astropliysics is to understand the magnetism 
of the Sun, and in so doing to develop the tools needed for the exploration of the magnetic 
activity in other types of stars across the Hertzsprung-Russel diagram. A crucial step toward 
this goal is to decipher the magnetic structure of the upper chromosphere, the key interface 
region between the underlying (cooler) photosphere and the overlying (hotter) corona where 
dominance of the physics passes from hydrodynamic to magnetic forces (e.g., the review by 
Harvey 2006). To this end, we need: 

• to identify observables sensitive to the magnetic fields of the upper chromosphere and 
transition region; 

• to develop suitable diagnostic tools to infer the magnetic field from those observables; 

• to design and build the telescopes and instrumentation needed for measuring the ob- 
servables. 

The upper chromosphere is indeed a very important part of the Sun because all of the 
plasma, magnetic field and energy in the corona and solar wind are supplied through this 
key boundary region. The observables that contain information on the physical conditions 
of this region are mainly spectral lines in the ultraviolet (UV) and far-UV (FUV) spectral 
range, with Lya being the most prominent line. The intensity profile (i.e., the Stokes /(A) 
profile) of Lya and other lines of the Lyman series have been measured on the solar disk 
by several instruments on board rockets and space-based telescopes (e.g., Roussel-Dupre 
1982; Warren et al. 1998). These measurements show that when observed on the solar disk 
the hydrogen Lyman lines are always in emission, and that the emission originates in the 
upper chromosphere and transition region. Interestingly, high-resolution Lya intensity im- 
ages of quiet-Sun regions taken with the Very high Angular resolution ULtraviolet Telescope 
(VAULT) sounding rocket show a forest of elongated fibrils suggesting closed magnetic loops 
in the upper chromosphere (Vourlidas et al. 2010). However, no quantitative information 
on the magnetic fied of the upper solar atmosphere can be obtained from this type of ob- 
servations because the intensity profile of Lya and of many other transition region lines is 
practically insensitive to the strength, inclination and azimuth of the magnetic field vector. 

The only way to obtain quantitative empirical information on the strength and orienta- 
tion of the magnetic field is through the measurement and interpretation of the polarization 
that some physical mechanisms introduce in spectral lines. Unfortunately, the familiar Zee- 
man effect is of little practical interest for the "measurement" of magnetic fields in the upper 
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chroinosphere and transition region (except perhaps in sunspots) because the circular polar- 
ization (Stokes V) scales with the ratio, 7^, between the Zeeman splitting and the Doppler 
line-width and the linear polarization (Stokes Q and U) scales with TZ^; this ratio is very 
small for the UV and FUV lines that originate in the weakly magnetized plasma of the upper 
solar chromosphere because TZ oc \B jyT (with A the wavelength, B the magnetic strength 
and T the kinetic temperature; see Landi Degl'Innocenti & Landolfi 2004). Fortunately, 
there is yet another physical mechanism by means of which the magnetic fields of a stellar 
atmosphere leave fingerprints in the spectral line polarization: the Hanle effect. The absorp- 
tion of anisotropic radiation produces atomic level alignment (i.e., the individual magnetic 
substates of levels with angular momentum j>l are unevenly populated, in such a way that 
the populations of sublevels with different values of \M\ are different and, moreover, quan- 
tum coherences between them may appear). This, in turn, gives rise to linear polarization 
in the spectral line under consideration (e.g., Trujillo Bueno 2001). The Hanle effect due to 
the presence of a magnetic field inclined with respect to the symmetry axis of the pumping 
radiation field modifies the atomic level alignment and the linear polarization of the spec- 
tral line radiation. Approximately, the emergent linear polarization is sensitive to magnetic 
strengths between 0.2Bh and 5Bh, where Bh= 1-157x10^^ /tincg is the critical Hanle field 
(in G) for which the Zeeman splitting of the line's level under consideration is similar to its 
natural width (with tufc the level's radiative lifetime, in seconds, and g its Lande factor). 
Note that Bh^50 G for Lya, Bh^^O G for Ly/3 and Bh^8 G for Ly7, and that the magnetic 
sensitivity provided by the Hanle effect is independent of the Doppler line width. 

Each Lyman line results from two blended transitions between a lower level, with j = 
1/2, and two upper levels, with j = 1/2 and j = 3/2. Therefore, for each Lyman line the 
only level that can be aligned and contribute to the emergent linear polarization is the upper 
level with j = 3/2 (when the hyperfine structure of hydrogen is neglected, which is known to 
be a good approximation for Lya). A particularly important question is whether the atomic 
alignment of the 2p^P3/2 level of Lya is significant or, in other words, if the Lya radiation that 
illuminates the hydrogen atoms of the solar atmosphere is sufficiently anisotropic in the line 
formation region. Obviously, the anisotropy of the incident Lya radiation is very significant 
in the optically thin case of the solar corona observed off-the-limb at large distances above 
it (e.g., Trujillo Bueno et al. 2005; and more references therein). The key question for 
us here is however whether the Lya radiation is sufficiently anisotropic within the solar 
transition region itself, so as to produce there a significant amount of atomic alignment in 
the 2p^P3/2 level. As mentioned above, the Lyman radiation that we observe when pointing 
to the solar disk is produced by the plasma of the solar transition region, which is very 
optically thick at the central wavelength of Lya. As a matter of fact, the center-to-limb 
variation (CLV) of the observed Lya intensity is very small or negligible (e.g., Roussel- 
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Dupre 1982; Curdt et al. 2008), which imphes that the intensity of the outgoing radiation 
(i.e., that propagating outwards) shows no significant variation with the hehocentric angle 6. 
However, this does not mean that the illumination of the transition region atoms is isotropic 
- which would imply that the scattering polarization in Lya is zero. As we shall see below, 
the intensity of the incoming radiation (i.e., that propagating inwards) calculated in semi- 
empirical and hydrodynamical models of the solar chromosphere shows significant CLV, so 
that the hydrogen atoms of the solar transition region must be illuminated by an anisotropic 
radiation field. The main goal of this Letter is to show that the linear-polarization amplitude 
of the emergent Lya radiation is predicted to be significant because of the anisotropy of the 
incoming radiation and that via the Hanle effect it is sensitive to the magnetic fields expected 
for the upper chromosphere and transition region of the Sun. 

We focus here on Lya because it is the most intense line of the solar transition region 
and also because its polarization is practically insensitive to collisional depolarization there, 
but we point out that linear polarization due to atomic level alignment is also expected for 
at least Ly/3 and Ly7. 



2. Formulation of the Problem 

We have applied the quantum theory of spectral line formation described in Chapter 7 
of Landi Degl'Innocenti & Landolfi (2004), which treats the scattering line polarization phe- 
nomenon as the temporal successsion of Ist-order absorption and re-emission processes, in- 
terpreted as statistically independent events. This complete frequency redistribution (CRD) 
approximation is suitable for reaching the objective of this Letter - that is, to estimate the 
line-center polarization amplitudes in the absence and in the presence of the Hanle effect 
(e.g., Sampoorna et al. 2010, and more references therein). 

The splitting between the 2p^Pi/2 and 2j9^P3/2 upper levels of Lya is of the order of 
10^° s~^, which is much larger than the level's natural width (of the order of 10^ s~^) but 
much smaller than the line's Doppler width (of the order of 10^^ s~^)- This justifies to neglect 
quantum interferences between such two j-levels (see Belluzzi & Trujillo Bueno 2011). The 
upper level 2j9^P3/2, with Lande factor g = 4/3, is the only level that contributes to the linear 
polarization of Lya when its hyperfine splitting is neglected, which as shown by Bommier & 
Sahal-Brechot (1982) is a good approximation. 

In order to estimate the linear polarization amplitudes of Lya we need to take into 
account radiative transfer effects in models of the solar atmosphere. To this end, we have 
chosen both the quiet-Sun semi-empirical model of Fontenla et al. (1993; hereafter the FAL- 
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C model) and the chromospheric hydrodynamical model of Carlsson & Stein (1997; hereafter 
the hydrodynamical model). The atomic model we have used includes all the fine-structure 
levels of the first three n-levels of hydrogen (that is, it includes Ly/3 and Ha, in addition to 
Lya). We quantify the excitation state of each j-level by means of the multipolar components 
of the atomic density matrix, whose self-consistent values at each spatial grid point of the 
model atmosphere have to be obtained by solving jointly the statistical equilibrium equations 
for such density-matrix components and the Stokes- vector transfer equation for each of the 
allowed radiative transitions in the atomic model. To this end, we have applied an efficient 
and accurate radiative transfer method (see Appendix A of Stepan & Trujillo Bueno 2011). 
Isotropic collisions with protons and electrons have been also taken into account, as described 
in Appendix B of Stepan & Trujillo Bueno (2011). In particular, the dipolar collisional rates 
between nearby j'-levels (i.e., between levels nlj and nl ± Ij ) have been calculated applying 
the semi-classical impact approximation theory (e.g., Sahal-Brechot et al. 1996), while the 
inelastic collisional rates between the different n-levels were taken from Przybilla & Butler 
(2004). We have included also the Stark broadening of the hydrogen lines, following Stehle 
(1996). With these physical ingredients, our non-LTE syntheses of the hydrogen intensity 
profiles are in good agreement with those computed by other researchers using more n-levels. 



3. The anisotropy of the Lya radiation in the solar transition region 

In the absence of magnetic fields the line-center amplitude of the fractional scattering 
polarization in the Lyman lines can be estimated by applying the following approximate 
expression (see Trujillo Bueno & Manso Sainz 1999, and note that the 1/v^ factor here 
takes into account that two blended transitions contribute to the intensity profile): 

7^^(i-/^')W2(j^j.)|, (1) 

where W2{ji,ju) = 1/2 for the Lyman lines, fi = cos6 (with 6 the heliocentric angle), and 
Jq/Jq is the fractional anisotropy of the spectral line radiation calculated at the height in 
the model atmosphere where the line center optical depth is unity along the line of sight 
(LOS). Therefore, the key quantity that determines the scattering polarization signals we 
may expect for the Lya line when observing the solar disk is Jq/Jq, where 

J'q = Jdu^^<p{u,n)i,^ (2) 

is the familar frequency-integrated mean intensity (0(z^, fl) being the normalized absorption 
profile, with u the frequency and O the ray direction), and 
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quantifies the anisotropy of tlie spectral line radiation tliat illuminates each point of the 
astrophysical plasma under consideratioiu. In solar-like atmospheres the first term of Eq. (3) 
plays the dominant role, so Jq>0 if the incident radiation field is predominantly vertical while 
Jq<0 if it is predominantly horizontal. The specific intensity depends on the spatial variation 
of the corresponding source function component, Sj. If Sj decreases with height in the 
solar atmosphere then the outgoing radiation shows limb darkening (i.e., it is predominantly 
vertical), while the incoming radiation shows limb brightening (i.e., it is predominantly 
horizontal). It is the magnitude of the vertical gradient of Si what determines the sign of 
Jg (see Fig. 4 of Trujillo Bueno 2001). For a constant source function we have that Jq<0, 
because in this particular case all outgoing rays have the same intensity and we are left only 
with the (predominantly horizontal) contribution of the incoming radiation. 

The left panel of Fig. 1 shows the spatial variation of Jq/Jq for Ha, Lya and Ly/3, 
calculated in the FAL-C model. Note that the anisotropies of the Lyman lines are practically 
zero all through the model atmosphere, except in the model's transition region where they are 
negative and significant (i.e., of the order of a few percent at the atmospheric heights where 
the line-center optical depth is unity along the LOS). The right panel of Fig. 1 illustrates 
that in the model's transition region the Lya anisotropy is dominated by the CLV of the 
incoming radiation. 



4. Scattering polarization and Hanle effect in Lya 

The two panels of Fig. 2 show the calculated center-to-limb variation of the emergent 
Q{\)/I{X) profile of the Lya line (hereafter, Q/I). The left panel corresponds to the results 
of our radiative transfer calculations in the semi-empirical FAL-C model, while the right 
panel shows the Q/I profiles that result from averaging the time- dependent Stokes I and Q 
profiles calculated in the hydrodynamical simulation model of Carlsson & Stein (1997). The 
variation with fi of the line center amplitude of the Q/I profile can be easily understood 
by using Fig. 1 (right panel) and Eq. (1). Remarkably, the amplitude and shape of the 
Q/I profiles calculated in both solar atmospheric models are qualitatively similar; this re- 
inforces our conclusion that scattering processes in the solar transition region should indeed 



^Note that these are the CRD expressions for J^ (with K = 0,2). Partial frequency redistribution 
expressions for J^ can be found in the paper by Sampoorna et aL (2010), which clarifies that CRD is a 
suitable approximation for estimating the line-center scattering polarization. 
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produce measurable linear polarization signals in Lja. We point out that the transition 
region starts at a rather precise height in the FAL-C model (i.e., at about 2200 km), while 
in the hydrodynamical model it fluctuates with time. 

We turn now to estimating the sensitivity of the emergent Lya linear polarization to 
the presence of organized and random magnetic fields in the solar transition region. Since 
similar conclusions are obtained with the hydrodynamical model, it suffices with showing the 
results of our self-consistent radiative transfer calculations in the FAL-C model atmosphere. 

Figure 3 shows the calculated Q/I and U/I signals for a LOS with /x = 0.3, which 
implies a distance of about 45" from the solar limb. Except for the B = OG case (shown in 
Fig. 3 with a gray dotted line for reference), the curves show Q/I and U/I for a horizontal 
(i.e., parallel to the solar surface) magnetic field of 20 G; the field has a fixed orientation 
given by the azimuth angle written in the figure panels. Note that both Q/I and U/I are 
sensitive to the azimuth of the magnetic field vector: while the Q/I signal changes but 
remains always negative, the sign of U/I also depends on the azimuth value. 

In such close-to-the-limb scattering geometry, U/I would be zero if, instead of having 
a dominant orientation, the magnetic field azimuth was uniformly distributed within the 
spatio-temporal resolution element of the observation. However, away from the solar disk 
center, Q/I remains sensitive to the magnetic field strength even in such an unfavorable sit- 
uation. This can be seen clearly in the left panel of Fig. 4, which shows how the amplitude 
of Q/I signal decreases as the magnetic strength of a random-azimuth horizontal field in- 
creases. As seen in this figure, for the scattering geometry of a close to the limb observation 
(e.g., fi = 0.3), the Hanle effect depolarizes. However, for the forward-scattering geometry of 
a disk center observation (i.e., a LOS with /i = 1), the Hanle effect of an inclined magnetic 
field with a given orientation creates linear polarization (see the right panel of Fig. 4 and 
note that in order to be able to detect the presence of a horizontal magnetic field of 30 G in 
the bulk of the solar transition region the polarimetric sensitivity of the measurement should 
be at least 0.1%). Therefore, the Hanle effect in Lya can be used to diagnose the magnetic 
field all over the solar disk (i.e., it is not restricted to an annular region around the solar 
limb). 



5. Conclusions 

Our radiative transfer investigation in semi-empirical and hydrodynamical models of 
the solar atmosphere indicate that the Lya line should show measurable linear polarization 
signals when observing the solar disk. Via the Hanle effect the line-center amplitudes turn 



out to be sensitive to the presence of magnetic fields in the solar transition region, with 
good sensitivity to field strengths between 10 and 100 G. For such field strengths there is no 
significant impact of the Zeeman effect on the Lya radiation. 

The above-mentioned prediction is based on the same quantum theory that explains the 
linear polarization profiles observed in various visible and near-IR lines (e.g., Manso Sainz 
& Trujillo Bueno 2003; Stepan & Trujillo Bueno 2010). This CRD theory is suitable also for 
estimating the line-center polarization amplitudes of strong resonance lines like Lya, but we 
point out and advance here that partial frequency redistribution (PRD) effects do produce 
interesting and sizable Q/I signals in the Lya wings (qualitatively similar to those shown in 
figure 9 of Trujillo Bueno 2011 for the Mg ll k-line at 2795 A). 

The predicted line-center amplitudes of the Lya scattering polarization are smaller than 
1%. Nevertheless, with a moderate aperture telescope they can be measured with a polari- 
metric sensitivity of 0.1% and a spectral resolution of 0.1 A by opting for a spatial resolution 
of the order of 1 arcsecond and a temporal resolution of the order of 1 minute (e.g., see 
the Lya polarization sounding rocket project outlined in Ishikawa et al. 2011). Informa- 
tion on the thermal and magnetic structure of the solar transition region can be obtained 
from the observed Q/I and U /I profiles themselves and through detailed radiative trans- 
fer simulations in given atmospheric models. Therefore, it would be worthwhile to develop 
a spectropolarimeter for measuring from space the linear polarization profiles caused by 
scattering processes and the Hanle effect in the Lya line of the solar transition region. Un- 
fortunately, the pioneering Soviet satellite experiment of Stenflo et al. (1980) failed because 
the optical transmission of the instrument was severely degraded due to in-orbit contamina- 
tion from the satellite. Such a development would open up a novel diagnostic window for 
the exploration of the magnetism of the outer solar atmosphere. 
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Fig. 1. — Tlie fractional anisotropy of the spectral line radiation. Left panel: the variation 
with height in the FAL-C model atmosphere of the CRD fractional anisotropy of the Lya 
(solid line), Ly/? (dashed line) and Ha (dotted line) radiation, with the arrows indicating the 
heights where the line center optical depth of each Lyman line is unity along line of sights 
with /i = 1 (left arrow) and yU = 0.1 (right arrow). Right panel: like in the left panel, but 
only for Lya and showing also the contributions of the outgoing radiation (with 0</i<l) 
and incoming radiation (with — l</i<0). 
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Fig. 2. — The unmagnetized case. Left panel: CLV of the emergent Lya Q/I profile calcu- 
lated in the FAL-C model. Right panel: CLV of the emergent Lya Q/I profile that results 
after time averaging the Stokes / and Q profiles computed at each time step of the hydro- 
dynamical model. The positive reference direction for Stokes Q is the parallel to the nearest 
limb. 
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Fig. 3. — The magnetized case. The Lya Q/ 1 and U /I signals for a LOS with /x = 0.3 
produced by scattering processes in the transition region of the FAL-C model atmosphere, 
taking into account the Hanle effect of a 20 G horizontal magnetic field. The various profiles 
correspond to the indicated values of the magnetic field azimuth (xb), measured counter- 
clockwise with respect to the projection of the LOS onto the solar surface plane. The positive 
reference direction for Stokes Q is the parallel to the nearest limb. 
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Fig. 4. — The magnetized case. In both panels the various Q/I profiles correspond to 
the indicated values of the magnetic field strength. Left: The Lya scattering polarization 
signals for a LOS with fx = 0.3 assuming the presence in the FAL-C model of a horizontal 
magnetic field with a random azimuth within the observational resolution element. The 
positive reference direction for Stokes Q is the parallel to the nearest limb. Note that in 
this close to the limb scattering geometry the linear polarization amplitude decreases as the 
magnetic strength is increased. Right: The Lya scattering polarization signals for a LOS 
with n = 1 assuming the presence of a horizontal magnetic field with a fixed azimuth. The 
positive reference direction for Stokes Q is along the magnetic field. Note that in this forward- 
scattering geometry the linear polarization amplitude increases as the magnetic strength is 
increased. 



